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ABSTRACT 

We look at the impact of Infra-Red (IR) and radio point sources on upcoming large yield Sunyaev- 
Zel'dovich (SZ) cluster surveys such as APEX, SPT and ACT. The IR and radio point source counts 
are based on observations by the the SCUBA and WMAP instruments respectively. We show that the 
contributions from IR source counts, when extrapolated from the SCUBA frequency of 350 GHz to the 
operating frequencies of these surveys (~ 100 — 300 GHz), can be a significant source of additional 'noise' 
which needs to be accounted for in order to extract the optimal science from these surveys. Alternatively, 
these surveys give us an opportunity to study IR sources, their numbers and clustering properties, opening 
a new window to the high rcdshift universe. For the radio point sources, the contribution depends on 
a more uncertain extrapolation from 40 GHz but is comparable to the IR near 200 GHz. However, the 
radio signal may be correlated with clusters of galaxies and have a disproportionately larger effect. 

Subject headings: cosmic microwave background — galaxies: clusters — cosmology: theory 



1. INTRODUCTION 

The study of anisotropies in the Cosmic Microwave Back- 
ground (CMB) has proven to be a gold-mine for cosmology. 
The primary anisotropies on scales larger than 10' have 
now been probed with high fidelity by WMAP (Bennett et 
al. (2003)) over the whole sky, leading to strong constraints 
on our cosmological model. Within the next few years this 
activity will be complemented by high angular resolution, 
high sensitivity observations of secondary anisotropies by 
the SZA\ APEX-SZ experiment^ the South Pole Tele- 
scope (SPT^) and the Atacama Cosmology Telescope (ACT^) 
which are aiming to make arcminute resolution maps with 
10/iK sensitivity (or better) at millimeter wavelengths. 

The dominant secondary anisotropy is expected to be 
the Compton scattering of cold CMB photons from hot 
gas along the line of sight, known as the thermal Sunyaev- 
Zel'dovich (SZ) effect (Sunyaev & Zel'dovich (1972; 1980); 
for recent reviews see Rephaeli (1995) and Birkinshaw 
(1999)) though other signals can be present at lower am- 
plitudes. In principle, a measurement of such anisotropies 
would constrain cosmological parameters (Weller et al. (2001), 
Levin et al. (2002), Majumdar & Mohr (2003), Hu (2003)), 
probe the thermal history of the intra-cluster medium (Ma- 
jumdar (2001), Zhang & Pen (2002)), put useful constraints 
on re-ionization models through the kinetic SZ effect (Zhang 
et al. (2003)) and allow us to map the dark matter back to 
the surface of last scattering (Seljak & Zaldarriaga (1999); 
Zaldarriaga & Seljak (1999); Hu (2001); Hirata & Seljak 
(2003); Okamoto & Hu (2003)). 

However, like the primary anisotropies, the secondary 
signals must be disentangled from other astrophysical emis- 
sions at the observed frequencies. One of the major source 
of contamination on small scales are point sources (by 
which we shall mean throughout sources unresolved at ar- 
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cminute resolution). The purpose of this paper is to out- 
line how point sources impact high resolution CMB exper- 
iments and the associated 'confusion' noise, describe how 
one can make small extrapolations of existing data to es- 
timate the amplitude of this noise and finally discuss the 
science which can come from studies of this source popula- 
tion. We make use of new observational constraints on the 
number of sources at frequencies and flux level relevant to 
APEX-SZ, SPT and ACT. 

2. POINT SOURCES AS NOISE 

We will see that point sources will be particular trou- 
blesome for high resolution surveys, and we will attempt 
to quantify their effect in terms of an equivalent "noise" . 
Such a quantification can only be approximate, however 
it is a useful metric by which to judge the relative impor- 
tance of different components and to plan survey strate- 
gies. Physically, the IR point sources are thought to be 
dusty, high redshift galaxies and so we do not expect them 
to be correlated with any particular place on the sky or any 
particular low redshift source (such as a cluster of galax- 
ies). We note, however, that the counts are sufficiently 
steep that lensing by the larger clusters can non-trivially 
increase the number of sources observed above a given flux 
cut and this should be taken into account in interpreting 
our results (see e.g. Perrotta et al. (2003) for recent theo- 
retical modelling). The radio sources are more problematic 
(Holder (2002)), as they are likely to be correlated with the 
SZ signal from clusters of galaxies to some extent. Since 
this correlation is difficult to quantify at present we will 
simply provide "noise" estimates for them as well. 

To proceed, we calculate the angular power spectrum 
for a population of point sources described by a flux dis- 
tribution and a clustering amplitude. If we assume that 
the number of sources of a given flux is independent of 
the number at a different flux, and if the angular two- 
point function of the point-sources is w(0), then the an- 
gular power spectrum, Cg, contributed by these sources is 
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(Scott & White (1999)) 

/•Scut J AT 

CKi^) = y^ 52 — d5.+«;,(^.f, (1) 

assuming that all sources with S>Scut are removed. Here, 
I,y = J S dN/dS dS is the background contributed by sources 
below S'cut- Following the conventional notation, Ci is the 
Legcndre transform of the correlation function C{9) pro- 
duced by the sources and we is the Legendre transform of 

C{e) = ^ 1^(2^ + l)CePe{cos 0) (2) 

e 

w{e) = ^Y.(2£+l)wePe{cose), (3) 

i 

with Pe{cos9) the Legcndre polynomial of order i. The 
first term in equation (1) is the usual Poisson shot-noise 
term (sec Peebles (1980) §46, or Tegmark & Efstathiou 
(1996)), the second is due to clustering, assuming that the 
clustering is independent of flux. All reasonable dN/dS 
give a Ci which convergers at the faint end and is very 
insensitive to the upper flux density cut. 

Because these sources are contributing a foreground to 
CMB anisotropy experiments we will recast our results in 
'temperature' units. This is easily done by applying the 
conversion factor between flux and temperature which is 
given by 
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where B^, is the Planck function, k is Boltzmann's con- 
stant, X = hi>/kBTcMB — i^/56.84GHz is the 'dimen- 
sionless frequency' and IJy = 10"^^ W m~2Hz~^. Conve- 
niently a 10 /iK CMB fluctuation in a 1' FWHM pixel gives 
a flux close to ImJy viz: 
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indicating that these upcoming surveys will be probing the 
source population at the mJy level. 

Since the clustering of the sources is highly uncertain^ 
at present, we will set = 0, though we shall return to 
this issue in §5. For now this is a conservative assumption 
as we expect the soiirces to be non-trivially clustered if 
they are associated with rare, highly biased tracers of the 
density field at high- 2; but we also expect those correlations 
to be most important at low-£. In the Poisson limit the 
Ci arc independent of I, i.e. they have the same shape as 
Poisson/white noise. We therefore express our results in 
terms of an 'eff'ective' point source noise CTpix per pixel of 
FWHM 6»pix using 

Ct = (^pixfpix)^ , (8) 

® Though we expect this to improve soon with the release of data 
from the SHADES survey, http://www.roe.ac.uk/ifa/shades/, and 
of course APEX-SZ, SPT and ACT. 



and express dpix in /kK for a fiducial 1' pixel. This value 
can be rescaled to any desired resolution using Eq. 8. This 
CTpix is closely related to the standard 'confusion noise' of- 
ten quoted by radio astronomers. To make the connection 
more explicit let us write 
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which defines g{x). Then the 'noise' induced by the sources 

is CTconf = Sq^N~,T{x) with X = Scut/So, Npi^ = A^O^'pix 

and 



ds s^(?(s) 
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If we assume we can subtract sources brighter than na this 
simplifies to 

CTconf = 'S'o - where ——= nWNp^ . (11) 
Vn/ -l-K^') 

Upon converting between flux and temperature units and 
accounting for the flux cut this expression is the same as 
Eq. 8. We mention finally one additional complication 
which we can include. If the sources have a frequency 
spectral index distribution J^{p) and we assume different 
sub-populations are independent then we can generalize 
the above to 

iVpix J d/3 n0)Sl{l3)I^ {xm (12) 

where Sq{(3) indicates the value of So obtained by ex- 
trapolating the fiducial So to the required frequency using 
Su oc . This tends to be a relatively small effect. Even 
a 30% gaussian uncertainty in /3 with (/3) = (see §3 be- 
low) increases the confusion noise at 5 mJy by only 10% 
over the /? = value. For this reason we shall work with 
uniform, constant (3 from now on. 

3. NOISE ESTIMATES 

The remaining step in computing the effective noise con- 
tributed by point sources is thus a model for the counts. 
While it is theoretically quite difficult to predict these 
functions, we are fortunate to have observations at close 
to the relevant frequencies and fiux levels. It is thus more 
robust to extrapolate the observations than to start from 
an ah initio theoretical model. 

3.1. Radio sources 



For the radio point sources, a fit to the Q-band data 
from WMAP (Bennett et al. (2003)) can be written as 

-2.7 

(13) 



dN 
dS^ 
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So 



where Nq = 1200 deg"^ = 4 x 10*' sr^^ with an uncer- 
tainty of around 30% and Sq ~ 1 mJy. In obtaining these 
numbers we need to extrapolate from very high flux levels 
to mJy levels, and this extrapolation is sensitive to the 
slope of the distribution. There is evidence that the slope 
flattens from -2.7 to -2.2 at lower fluxes. We shall choose 
-2.3 as a compromise, since S^^^dN/dS is roughly flat in 
Fig. 13 of Bennett et al. ((2003)). For So = ImJy this 
lowers No to SOdeg"^. The extrapolation to higher fre- 
quency is also somewhat uncertain. The primary emission 
mechanism in these sources is synchrotron emission, which 
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for a power-law electron spectrum would give a power-law 
Sv However synchrotron 'aging' and optical depth effects 
can give departures from a power-law spectrum, allowing 
a wide range of spectral shapes including spectra which 
peak at or above 100 GHz. Advection dominated accretion 
onto a super-massive black hole is also characterized by a 
strongly inverted spectrum peaking at millimeter wave- 
lengths, as are radio afterglows of gamma-ray bursts. In 
addition, many of the sources with significant flux at high 
frequency mc thought to be; young, compact sources which 
are likely to be highly time variable. 

Given these large uncertainties, we will attempt to bracket 
the reasonable range and extrapolate the fluxes to higher 
frequency assuming two different power-law spectral in- 
dices, S„ (X with j3 — —0.3 (Tegmark & Efstathiou 
(1996)) and /? = (close to the mean of the distribution 
in Trushkin (2003)). At the mJy level we have close to 
one radio source per 50 beams if the point sources pre- 
dominantly have a flat (or rising) spectrum. 

3.2. IR sources 

For the IR sources we can use observations at 350 GHz 
made with the Submillimeter Common-User Bolometer 
Array (SCUBA; Holland et al. (1999)) on the James Clerk 
Maxwell Telescope. SCUBA has been used to make several 
deep observations (Barger et al. (1998); Eales et al. (1998); 
Holland et al. (1998); Hughes et al. (1998); Smail et al. (1997)) 
from which we can extract source counts. In particular we 
shall use the recent work of Borys et al. (2003) who give a 
phenomenological fit to the counts at 350 GHz: 

s=t[(i)Ki)T' 

with No = 1.5 X lO^deg-^ = 4.9 x lO^sr-i and So ~ 
1.8 mJy. This model provides a reasonable fit to the ex- 
isting data near 1 mJy that does not over-produce the 
far-infrared background (FIB) light (Puget et al. (1996)). 
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Fig. 1. — The cumulative SCUBA source counts at 350 GHz as 
a function of flux, S, from Smail et al. (1997) (squares) and Borys 
et al. (1999; 2003) (circles). The fit of Borys et al. (2003) is shown 
as the solid line, and the extrapolations to 220 GHz and 150 GHz 
assuming oc v^'^ as dashed lines. Upcoming surveys will be 
sensitive to sources with 5 ~ 1 mJy. 



We estimate that the uncertainty in the normalization is 
roughly a factor of 2, due to the small sky area surveyed. 
We will extrapolate from 350 GHz to lower frequency us- 
ing oc v^^ with (3 = 2.5, close to the typical spectrum of 
a dusty star-bust galaxy at high-z. We will also examine 
the effect of steeping the frequency dependence to a more 
conservative value /3 = 3. 

In figure 1, we have plotted the observed IR source 
counts of (Borys et al. (1999), Smail et al. (1997)) at 
350 GHz along with our extrapolation to lower frequen- 
cies. Note that at 150 GHz there would be more than 100 
IR point sources above ImJy per deg^! This means tens 
to hundreds of thousands of IR sources (depending on the 
survey). These sources need to be carefully accounted for 
when considering secondary CMB science. 

4. RESULTS 

For 150 and 220 GHz the point source contributions, as- 
suming 1' pixels, are summarized in Tables 1 and 2. We 
do not include the effects of source clustering, which would 
increase these numbers, or gravitational lensing of the pop- 
ulation by a foreground object such as a cluster. The corre- 
sponding power spectra are shown in Figs. 2 (at 150 GHz) 
and 3 (at 220 GHz) . We caution the reader that several 
of the processes on these figures arc non-Gaussian, so the 
power spectra tcill only some of the story. Nonetheless 
we can see that point sources will be a large contribution 
to the signal at these frequencies, which will need to be 
properly accounted for in any analysis aimed at recover- 
ing the thermal SZ effect. Even more care will be required 
for processes with lower signal levels or methods which 
assume that the underlying map is predominantly Gaus- 
sian in nature (as for example the primary CMB is) since 
residual emission needs to be carefully controlled. 

There is one obvious strategy for controlling the effect of 




Fig. 2. — The angular power spectra of various sources at 150 GHz. 
The line labeled "primary" is the primary CMB anisotropy spec- 
trum, the line labeled "tSZ" is the thermal SZ spectrum from White, 
Hernquist & Springel (2002) and is uncertain at the factor of 2 level. 
The (dashed) line labeled "noise" is the instrument noise assuming 
10 /iK per 1' pixel and a resolution of 1' . The two unlabeled lines 
rising rapidly to the top right of the plot are the Poisson contribu- 
tion of the radio and IR sources for 5cut = 5 mJy, assuming /3 = 
and 2.5 respectively. 
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Fig. 3. — The angular power spectra of various sources at 220 GHz, 
the thermal SZ "null" . Lines are as in Fig. 2 except that a model of 
the correlations in the FIB from Knox et al. (2001) is also shown. 
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Table 1 

The effective noise for Poisson distributed radio 

SOURCES AT 150 and 220 GHZ ASSUMING A FIDUCIAL l' 
PIXEL. We show EXTRAPOLATIONS FROM 40 GHZ ASSUMING 
A SPECTRAL INDEX OF /? = AND —0.3. 



these point sources, and that is to use multiple observing 
frequencies and multiple instruments with matched sensi- 
tivities and resolutions. 

Finally we make an estimate of the 'confusion noise' as 
a function of beam size at these frequencies. For the radio 
sources I(x) = -J x'^~'^-^ j — 2.3) allowing us to estimate 
CTconf as a function of beam size. If we assume /3 = the 
typical values for a 4cr cut at 150 GHz arc (Tconf = 2mJy 
at 5', 0.4 mJy at 2', 0.1 mJy at 1' and 0.05 mJy at 0.5'. For 
the IR sources T{x) is a hypergeometric function. As a; — > 
0, x/T{x) \/2 while x/T{x) oc a; for x 1. Assuming 
/3 = 2.5 the typical values for a 4it cut at 150 GHz are 
fconf = 3.8 mJy at 5', 1.4 mJy at 2', 0.6 mJy at 1' and 
0.3 mJy at 0.5'. In thermodynamic units this is e.g. 20 /uK 
at 1' in agreement with Table 2. We show the frequency 
dependence of the confusion noise in Fig. 4. 

5. DOING POINT SOURCE STUDIES WITH APEX/SPT 

The point sources described above represent an analysis 
challenge for studies of secondary CMB anisotropics such 
as the thermal and kinetic SZ effects or gravitational lens- 
ing. However, they also represent increased science reach 
for the surveys in different fields of research. 

At present, the SCUBA sources account for 40% of the 
350 GHz sub-mm background (Borys et al. (2003)) and 
thus go a long way towards resolving the high- 2; part of 
the cosmic far-infrared background (FIB). This suggests 
that APEX-SZ, SPT and ACT with their high sensitivity 
and large arcal coverage could become ideal instruments 
for studying the FIB at the faint but especially the bright 
end of the source counts. As is evident from Figure 1, these 
future surveys would also be able to detect thousands of 
high redshift sources. 

These upcoming surveys would be able to probe the 
FIB correlations which has been shown (Haiman & Knox 
(2000), Knox et al. (2001)) to give an additional handle on 
early structure formation at z > 1. As wc show in Fig. 3, 
detailed observations of the FIB correlations are possible if 
we can accurately subtract the contribution from primary 
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Table 2 

The effective noise for Poisson distributed IR 

SOURCES AT 150 AND 220 GHz ASSUMING A FIDUCIAL 1' 

PIXEL. We show EXTRAPOLATIONS FROM 350 GHZ 
ASSUMING A SPECTRAL INDEX OF /3 = 2.5 AND 3. AnY 
CORRELATIONS IN THE SOURCES WILL (LIKELY) INCREASE 
f^pix- 
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Fig. 4. — Confusion noise, in /nK, for a 1' beam as a function 
of frequency. Contributions from Poisson distributed radio and IR 
sources are shown. The solid lines represent /? = and 2.5 while the 
dotted lines represent P = —0.3 and 3 respectively. 
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anisotropics and measure power on angular scales of a few 
degrees. Under assumptions similar to those in Knox et 
al. (2001), we expect that a 10^ square degree survey to 
around 20 /xK would make precision measurements (at the 
several percent level near I ^ 10'^) of the clustering of the 
FIB sources. Multiple frequency information will be cru- 
cial in separating the signal from other contaminants. In 
addition, since the FIB is composed of contributions from 
sources at different redshifts, the FIB sky maps at differ- 
ent frequencies are not perfectly correlated. The shape of 
the FIB power spectrum at different frequencies and the 
correlations between them can give information about the 
contributing sources (Knox et al. (2001)). Studies of the 
FIB with these surveys would be complimentary to the 
high angular resolution optical and UV observation of the 
individual sources. 

For the radio sources, these surveys operate at higher 
frequencies and are complementary to low frequency com- 
pilations of radio point sources (Condon et al. (1998), 
White et al. (1997)), giving us the opportunity to study 
the spectral dependence of known sources and to probe 
an entirely new population. Since, in general, radio point 
sources have non-trivial spectra (Herbig & Readhead (1992)), 
the multi-frequency detection of such sources would be 
invaluable in our understanding of the behavior of radio 
sources. Under reasonable assumptions associating differ- 
ent populations of radio sources (of a given lifetime) with 
halos (Haiman & Hui (2001), Martini & Weinberg (2001)), 
one would be able to constrain cosmological models through 
a study of the clustering of radio sources. At 1 mJy we ex- 
pect about 1 source per 50 arcmin^, which would allow 
APEX-SZ to constrain l + w{d) on arcminute scales at the 
several times 10~^ level. Whether this is enough to per- 
form a reliable measurement oiw{9) depends on the degree 
of bias and line-of-sight dilution of the clustering. At lower 
frequencies, the angular correlation function w{0) oc 0^~'' 
with 7 ~ 1.8. The amplitude. A, of 'w{9) depends on the 
flux limit and flattens out for low flux limits. Extrapo- 
lating from the NVSS data one would expect A< 10~^ at 
1 mJy (Overzier et al. (2003)), which would be near or 
below the threshold for APEX-SZ. 

Clearly to extract all of the excellent science that can 
be done with these instruments will require coordinated 
observations over a range of frequencies. In this regard the 
sub-mm instruments intended for the APEX platform or 
the planned SCUBA-2'' with its higher angular resolution 
and frequency range will be particularly valuable. For the 
radio sources coordinated observations with the SZA or 
other lower frequency instruments would be useful. 

6. DISCUSSION AND CONCLUSIONS 

Advances in detector technology have brought us to the 
stage where sensitive, high angular resolution, wide area 
surveys in the sub-mm are practical. Several new instru- 
ments are funded and under construction which should 
map hundreds or thousands of square degrees of sky at ar- 
cminute resolution with sensitivities of around 10 /iK. The 
stated aim of these surveys is to find clusters of galaxies 
using the Sunyaev-Zel'dovich effect. 

At these sensitivities, foregrounds are of particular im- 
portance, and in this paper we have used new observa- 

^ http://www.roe.ac.uk/atc/projects/scuba_two/ 



tional data on the IR and radio source counts to esti- 
mate their impact on upcoming surveys. The contribution 
from both radio and IR point sources are found to be non- 
negligible. The IR point sources are of particular interest, 
because their number density is higher than was expected 
based on early theoretical modeling. Using a model of the 
source counts we predict that surveys with 1' resolution 
will be 'confusion limited' before they reach 10 /xK. Typ- 
ical effective noise levels from unsubtracted point sources 
can be more than twice the nominal survey sensitivities 
at 150 GHz and 1' resolution and even higher at 220 GHz. 
The contamination by point sources has a strong impact on 
how such surveys will study the thermal SZ effect, and may 
provide the ultimate limitation to studies of the kinetic SZ 
effect and weak gravitational lensing. Such surveys need 
to make multi-frequency observations or coordinate obser- 
vations with other instruments to allow them to remove 
point sources to low flux levels. 

Conversely these instruments will have a golden oppor- 
tunity to study the inverted spectrum radio sources un- 
derrepresented in current surveys, map the FIB and probe 
the high-2 universe. Not only will these upcoming surveys 
have the capability to detect thousands of high redshift 
IR and radio sources to very low flux limits, they will also 
be able to look into the correlations of the FIB and the 
clustering of radio point sources. 
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